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Preface 


This  report  documents  the  foundation  investigation  conducted  by  the  US 
Army  Engineer  Waterways  Experiment  Station  (WES)  for  the  Ground  Based  Radar 
Project  at  Kwajalein  Island,  US  Army  Kwajalein  Atoll,  Marshall  Islands.  The 
work  was  performed  during  the  period  25  October  through  3  November  1989  for 
the  Engineering  Division  of  the  US  Army  Engineer  Division,  Pacific  Ocean 
(POD),  under  IAO  No.  E87890083  POD  MIL-R  dated  8  Sep  89. 
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Laboratory  (GL) ,  WES,  was  the  Project  Engineer  for  this  study.  Mr  Michael  K. 
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FOUNDATION  INVESTIGATION  FOR  GROUND  BASED  RADAR  PROJECT 
KWAJALEIN  ISLAND,  MARSHALL  ISLANDS 

Introduction 

1.  Background.  The  foundation  investigation  was  performed  adjacent  to 
the  Defense  Control  Center  (DCC)  building  which  is  an  existing  structure  on 
which  a  tower  will  be  constructed  to  support  a  radar  antenna.  The  existing 
building  foundation  needed  to  be  evaluated  to  ensure  that  it  would  provide  a 
stable  foundation  during  the  dynamic  loading  of  this  new  antenna  support 
system.  Testing  within  the  building  itself  was  not  possible;  therefore,  tests 
were  conducted  near  the  building  in  an  area  that  was  included  in  the  earth 
berm  surcharged  foundation  of  the  existing  DCC  building.  It  is  assumed  that 
site  conditions  at  the  test  locations  and  beneath  the  building  are  virtually 
the  same  because  of  the  proximity  and  similar  preloading  conditions. 

2.  Purpose .  The  purpose  of  this  investigation  was  to  assess  the 
foundation  materials  at  the  site  to  a  depth  of  60  ft^  by  employing  in  situ 
geophysical  and  geotechnical  methods.  Geophysical  methods  were  used  to 
determine  compression  (P)-  and  shear  (S)-  wave  velocities  so  that  a  velocity 
zonation  of  the  foundation  materials  could  be  determined.  For  this  study,  a 

1.  A  table  of  factors  for  converting  non-SI  to  SI  (metric)  units  of 
measurement  is  presented  on  page  6. 
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suite  of  seismic  methods  consisting  of  surface  refraction,  downhole,  and 
crosshole  tests  were  conducted  to  determine  the  above  values.  The 
geotechnical  methods  employed  Standard  Penetration  Testing  (SPT),  density 
determinations,  gradation,  classification,  and  indexing  tests  of  selected  soil 
samples.  Also,  a  Plate  Bearing  Test  was  performed  at  one  location  at  the 
site.  The  geotechnical  tests  were  performed  by  POD  personnel.  The  SPT's  and 
laboratory  sample  testing  provided  N- values,  material  classif ication  and 
density  of  the  foundation  materials.  Knowing  this  information  the  elastic 
constants  Young's  modulus  (E) ,  Shear  modulus  (G) ,  and  Poisson's  ratio  (-J)  can 
be  determined  for  the  foundation  materials  which  are  needed  in  the  design  of 
the  antenna  support  structure.  The  plate  bearing  test  provided  an  alternate 
method  for  an  in  situ  determination  of  E. 

Site  Description 

3.  General .  The  location  of  this  study  was  the  US  Army  Kwajalein  Atoll 

(USAKA)  which  is  located  in  the  northern  Marshall  Islands  in  the  west  central 

Pacific  Ocean.  The  site  is  located  at  the  western  end  of  Kwajalein  Island 

adjacent  to  the  western  side  of  the  DCC  building,  as  shown  in  figure  1.  The 

regional  geology  of  the  site  is  an  island  derived  from  the  buildup  of  coral 

skeletons  on  the  submerged  rim  of  an  extinct  volcano.  The  topography  is  low 

o 

and  flat.  The  elevation  (el  )  of  the  site  is  approximately  8  ft  above  mean 

2.  All  elevations  are  in  feet  and  are  referenced  to  mean  sea  level 


8 


9 


Site  Location 


sea  level.  The  materials  composing  the  subsurface  are  unconsolidated 
calcareous  materials  on  a  reef  surface.  The  subsurface  materials  within  the 
landmass  are  composed  of  unconsolidated  limestone  derived  sediments  of  cobbles 
to  silt  sizes  with  a  general  USCS  classification  of  SP.  Also,  the  depth  of 
the  reef  rock  increases  rapidly  in  the  direction  of  the  lagoon  (US  Army 
Engineers,  1989).  This  general  geologic  description  of  the  site  agrees  well 
with  the  information  derived  from  the  borings  drilled  to  accommodate  the 
subsurface  geophysical  tests. 

4.  Foundation .  The  foundation  at  the  site  is  composed  of  naturally 
occurring  materials  below  a  depth  of  14  ft  (el  -6).  These  materials  are 
composed  of  sands  with  silt  and  gravel  with  the  occurrence  of  gravel 
increasing  with  depth.  The  boring  logs  contained  no  indication  of  coral 
limestone  rock  within  the  upper  60  ft.  This  information  was  obtained  from  the 
boring  logs  from  holes  BH-1,  BH-2,  and  BH-3  which  are  included  in  Appendix  A. 

A  fill  of  beach  sand  was  hydraulically  placed  to  a  height  of  8  ft  above  the 
original  island  surface.  This  fill  was  surcharged  with  8  ft  of  additional 
fill  which  was  removed  prior  to  start  of  construction.  The  water  table  is 
found  at  an  average  depth  of  8  ft  (el  0)  and  fluctuates  ^2  ft  with  the  tide. 

5.  Idealization.  For  clarity  and  use  in  engineering  analysis  it  is 
necessary  to  idealize  the  foundation  materials  into  discrete  layers  which  will 
then  be  identified  by  material  type.  Material  properties  are  assigned  to  each 
layer  based  upon  the  test  results.  The  geological  setting  and  construction 
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history  of  the  site  suggest  idealization  into  four  layers.  The  top  layer 
consisting  of  fill  will  be  variable  and  contain  a  shallow  near  surface  zone  of 
fill  that  may  be  disturbed  and  contain  various  construction  materials  buried 
during  the  construction  and  cleanup  of  the  site  in  conjunction  with  the 
construction  of  the  DCC  building.  This  zone  would  then  change  into 
undisturbed  hydraulic  fill  placed  above  the  water  table.  A  change  in 
character  of  the  fill  placed  below  the  water  table  might  be  evidenced  and 
would  be  caused  by  the  different  placement  conditions.  Layer  two  would  be  the 
naturally  occurring  near  surface  ocean  deposited  beach  materials.  The  third 
and  fourth  layers  would  consist  of  the  coarser  and  more  dense  unconsolidated 
materials . 


Test  Program 

6.  The  locations  of  tests  performed  during  this  investigation  are  shown 
in  figure  2.  All  phases  of  the  geophysical  test  program,  except  the  crosshole 
S-wave  test,  were  conducted  according  to  Engineering  Manual  1110-1-1802 
guidelines  (Department  of  the  Army,  1979).  The  surface  portion  of  the  test 
program  consisted  of  three  seismic  refraction  lines  (R1,R2,R3)  in  the  vicinity 
of  the  plate  bearing  test  location.  Crosshole  and  downhole  testing  were 
performed  near  the  DCC  building  in  three  boreholes  (BH-1,2,3)  spaced  10  ft  and 
5  ft  apart,  respectively,  which  were  drilled  to  a  depth  of  60  ft.  In  addition 
to  the  geophysical  tests  employed  in  the  completed  boreholes,  SPT  and  sampling 
were  performed  during  the  drilling  operation. 
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Figure  2.  Test  Layout 
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Geophysical  Tests  Procedures 

7.  Crosshole .  In  preparation  for  the  crosshole  testing,  8-in.  diameter 
boreholes  were  drilled  to  the  60- ft  depth.  These  holes  were  then  cased  with 
schedule  40  4- in.  ID  PVC  pipe  with  the  annular  space  between  the  borings  and 
casings  grouted  with  a  material  that  approximates  the  consistency  of  soil  when 
it  sets.  The  vertical  alignment  of  each  borehole  was  checked  to  see  if  there 
was  any  appreciable  drift  since  accurate  reduction  of  the  crosshole  data 
requires  that  the  straight-lined  distance  be  known  between  source  and  receiver 
at  each  test  elevation.  The  S-wave  test  procedure  consisted  of  placing  a 
downhole  triaxial  geophone  array  in  the  receiver  hole(s)  and  a  downhole 
electrically  powered  vibrator  in  the  source  hole.  The  vibrator  is  frequency 
and  duration  controlled  and  produces  a  repeatable  and  vertically  polarized 
shear  (SV)  -wave  which  allows  accurate  arrival  time  determination.  The  source 
and  receiver(s)  were  lowered  to  the  same  depth  in  a  borehole  set  (one  boring 
for  the  source  and  two  borings  for  the  receivers)  and  clamped  to  the  casing 
walls  using  inflatable  bladders  or  pneumatically  powered  rams.  The  vibrator 
frequency  was  varied  between  50  and  500  Hz  using  a  four  cycle  burst  mode  and 
monitored  until  an  optimal  frequency  was  found  that  propagated  well  at  that 
depth.  The  source  waveform  and  the  received  waveform  were  recorded  using  a 
digital  seismograph.  The  data  were  stacked  (enhanced)  until  a  well  defined 
waveform  was  produced.  For  the  P-wave  test,  the  seismic  source  was  an 
exploding  bridgewire  (EBW)  detonator  which  was  sufficiently  strong  in  energy 
that  data  stacking  was  not  necessary.  Several  different  test  configurations 
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were  employed  in  the  test  program  that  are  tabulated  in  figure  3.  Knowing  the 
distance  between  the  source  and  receiver  and  the  P-  and  S-  wave  arrival  times 
at  each  test  depth,  an  analysis  of  these  data  sets  was  made  with  the  aid  of  a 
computer  program  "CROSSHOLE"  (Butler  et  al,  1978).  This  program  calculates 
true  P-  and  S-wave  velocities  and  determines  velocity  zones  and  depths  to 
interfaces . 

8.  Downhole.  The  downhole  test  is  similar  to  the  crosshole  test  except 
the  source  is  kept  at  the  surface  while  the  receiver  array  is  lowered  in  a 
boring  at  5- ft  intervals.  The  source  for  the  S-wave  test  is  a  hammer  striking 
a  wooden  plank  on  alternate  ends,  which  produces  two  records.  The  seismic 
signals  produced  by  this  procedure  are  predominantly  horizontally 
polarized  S- waves,  with  polarity  depending  on  the  direction  of  the  hammer 
strike.  The  signals  detected  by  the  horizontal  geophones  on  these  two  records 
are  overlain  and  examined  for  a  polarity  reversal  which  is  considered  the 
arrival  of  the  S-wave.  The  P-wave  source  for  this  test  was  a  downward  hammer 
blow  to  a  steel  plate  with  the  vertical  geophone  signal  being  used  to 
determine  the  P-wave  arrival.  The  data  are  reduced  by  plotting  arrival  times 
versus  slant  distance  between  source  and  receiver.  The  inverse  slope  of  the 
line  segments  drawn  through  the  data  points  gives  the  velocities  and  slope 
changes  in  the  line  segments  indicate  the  approximate  depths  where  the 
velocities  change.  The  various  test  configurations  used  in  this  study  are 
shown  in  figure  3. 
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TEST  CONFIGURATIONS 


BOREHOLE  TESTS 


SURFACE  TESTING 

geophone  spacing,  ft  5,10,10,10,10,10,10,10,10,10,10,5 

;  R1  Ao  oooooooooo  CA 

-  110  ft  - 

j 

geophone  spacing,  ft  2, 2, 2, 2, 2, 2, 5, 5, 5, 5, 5, 5 

!  R2  Ao  OOOOOOOOO  GO 

i 

| - 42  ft  - 

geophone  spacing,  ft  2, 2, 2, 2, 2, 2, 3, 5, 5, 5, 5, 5 
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9.  Surface  seismic  refraction.  The  procedure  for  the  surface 
refraction  tests  was  to  place  geophones  along  a  straight  line  on  the  ground 
surface  at  selected  intervals  with  an  energy  source  initiated  at  a  selected 
offset  distance  from  each  end  of  the  geophone  array.  Two  types  of  surface 
tests  were  performed  which  consisted  of  a  P-wave  and  S-wave  test.  Lines  R1 
and  R2  were  P-wave  tests  and  line  R3  was  a  S-wave  test.  The  energy  sources 
for  these  tests  were  identical  to  the  ones  used  for  the  downhole  P-  and  S- 
wave  tests.  The  test  configuration  for  each  line  is  shown  in  figure  3.  The 
data  reduction  consisted  of  plotting  first  arrival  time  of  the  P-  and  S-  wave 
signal  detected  at  each  geophone  versus  geophone  distance  from  the  source. 
From  these  time  versus  distance  (TD)  plots,  velocities  and  depths  to 
refracting  interfaces  were  determined  using  the  computer  program  "SEISMO" 
(Yule  and  Sharp,  1989). 


Geotechnical  Test  Results 

10.  Soil  classification  and  density.  Analysis  of  the  soil  sampling 
data  has  produced  a  four  layer  profile  interpretation  which  is  illustrated  in 
figure  4.  The  first  layer  of  fill  material  is  a  medium  dense  to  dense  poorly 
sorted  sand  (SP)  with  zones  of  silty  sand  and  some  cobbles  and  gravel.  This 
layer  extends  to  a  depth  of  13  ft.  The  second  layer  begins  with  the  original 
soil  deposit.  This  layer  is  a  medium  dense  silty,  very  fine  grained  sand  (SM) 
which  can  be  found  between  14  and  19  ft  in  depth  (el  -6  to  -11).  The  third 
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FILL 


layer  is  a  medium  dense  well  graded  sand  with  some  silt  and  gravel  that  ends 
at  a  depth  of  29  ft  (el  -21).  The  fourth  layer  is  a  loose  to  medium  dense  to 
dense  silty  sand  with  an  increasing  gravel  content  below  50  ft  (el  -42) . 

11.  Density  determinations  of  samples  above  the  water  table  in  the 
fill  layer  show  an  average  moist  unit  weight  of  100  lb/ft^  and  a  water  content 
of  11  percent.  The  specific  gravity  of  the  solids  was  measured  to  be  2.8. 
These  results  were  used  to  estimate  a  total  unit  weight  of  the  soil  below  the 
water  table  of  120  lb/ft  .  Detailed  information  of  the  testing  program  and 
individual  sample  test  results  can  be  found  in  the  boring  logs,  Appendix  A. 

12.  SPT.  The  SPT's  performed  in  boreholes  BH- 1,2, and  3  are  presented 
in  figure  5  as  two  plots  of  N-values  versus  depth.  The  recorded  SPT  data  can 
be  found  in  the  boring  logs,  Appendix  A.  The  leftmost  plot  presents  the  raw 
or  measured  blowcounts.  In  the  next  plot,  the  (Nl)gg  value  corresponding  to 
each  measured  N-value  is  shown.  The  (Nl)gg  value  is  determined  by  adjusting 
the  observed  N-value  to  an  equivalent  N-value  had  the  test  been  conducted  at  a 
vertical  effective  overburden  stress  of  1  ton/ft  at  an  energy  level  of  60 
percent  of  the  theoretical  maximum  applied  energy  of  the  drop  weight.  The 
SPT's  from  this  project  were  performed  using  the  rope  and  cat-head  system 
which  operates  at  the  60  percent  energy  level.  Determining  an  equivalent 
<N1>60  blowcount  allows  the  results  to  be  compared  with  other  equivalent 
(Nl)go  values  and,  therefore,  established  correlations  of  relative  density  can 
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MEASURED  EQUIVALENT  BLOWCOUNTS 
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Figure  5.  Standard  Penetration  Test  Results 


was 


be  used  to  estimate  these  parameters  for  the  soil  at  the  site.  The  (Nl)gg 
obtained  by  multiplying  the  measured  blow  counts  by  the  correction  factor  Cn, 
which  was  determined  from  the  curves  presented  in  figure  6.  The  vertical 
effective  stress  was  calculated  using  a  total  unit  weight  of  100  lb/ft  for 

"3 

the  soil  above  the  water  table  and  a  saturated  unit  weight  of  120  lb/ft  for 
soil  below  the  water  table.  A  depth  to  the  water  table  of  8  ft  was  used  in 
these  calculations. 

13.  The  SPT  data  show  considerable  scatter  especially  near  the 
surface.  The  average  measured  and  average  equivalent  (N1)^q  blowcounts  are 
annotated  in  figure  5  and  in  table  1  for  each  soil  layer.  These 
representative  values  for  each  layer  are  plotted  as  vertical  lines  on  the 
plots.  Using  the  (Nl)gg  values  and  the  correlation  curve  for  sands  given  in 
figure  7,  estimates  of  the  relative  density  of  the  sands  at  the  site  were 
determined  and  listed  in  table  1.  A  site  global  average  (N1)^q  value  of  22  is 
appropriate  for  the  sands  and  correlates  to  sand  with  a  relative  density  of  65 
percent . 
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Relationship  Between  Relative  Density  and  (Nl)fe0  of  Sands 


Table  1.  SPT  Test  Results 
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1 

“T" 

N(measured) 

<N1>60 

(based  on  figure  7)  ! 

i 

i 

1  Tokimatsu  &  Seed  1987) 

I 

1 

20 

25 

67 

2 

17 

21 

64 

3 

I  21 

21 

64  ; 

4 

I  20 

15 

t  1 

55  ; 

14.  Plate  Bearing  Test,  A  Plate  Bearing  Test  was  conducted  at  one 
location  and  consisted  of  static  and  cyclic  *_  ...  tests  on  a  26.6-  in.  square 
steel  plate.  The  testing  elevat'on  was  1  ft  below  the  existing  grade  (el  7). 
The  complete  description  of  these  tests  and  their  results  are  presented  in 
Appendix  B.  For  the  static  plate  bearing  test,  a  maximum  load  of  34.9  lb/in. 
was  applied.  Based  upon  this  test,  a  modulus  of  soil  reaction  of  340  lb/in. 

was  determined.  The  cyclic  plate  bearing  test  was  performed  with  a  19.3 

2  2 
lb/in.  static  load  and  a  cyclic  load  of  5  lb/in.  .  From  this  test  an  E  of 

3  2 

123.6x(10  )  lb/in.  was  determined.  This  computation  is  presented  in 
Appendix  B. 
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Geophysical  Test  Results 

Surface  Seismic  Refraction 

15.  P-wave  Tests  R1  and  R2 .  The  results  from  P-wave  surface 
refraction  lines  R1  and  R2  are  shown  in  figure  8.  The  R2  line  was  configured 
to  investigate  the  near-surface  materials;  whereas,  the  longer  R1  line  was 
employed  to  investigate  to  deeper  strata.  The  two  overlapping  data  sets  were 
combined  to  provide  a  detailed  interpretation.  From  this  composite  TD  plot  a 
four  layer  seismic  profile  was  indicated.  A  cross  section  is  also  shown  in 
figure  8  with  the  depths  to  interfaces  and  velocities  presented.  The 
calculated  depths  at  each  end  of  the  line  indicate  very  little  dipping  of  the 
layers  in  this  av  i.  Layers  one  and  two  with  velocities  of  1,440  ft /sec  and 
2,610  ft/sec,  respectively,  correspond  to  the  fill  materials  above  the  water 
table.  The  third  layer  with  a  velocity  of  5,110  ft/sec  is  the  saturated  fill 
material.  This  result  shows  the  "seismic"  water  table  to  be  at  an  average 
depth  of  8.0  ft  (el  0).  The  fourth  layer  detected  at  a  depth  of  17-19  ft 
(el  -9  to  -11)  has  a  velocity  of  6,400  ft/sec  and  corresponds  to  the  saturated 
well  sorted  sands  approximately  6  ft  below  the  bottom  of  the  fill. 

16.  S-wave  Test  R3 .  Results  from  the  surface  shear  wave  refraction 
line  R3  are  not  presented  because  the  P-wave  arrivals  were  so  strong  that  the 
S-wave  arrivals  could  not  be  accurately  identified  and  analyzed. 
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Surface  Seismic  Refraction  Test  Results  Lines  Rl  atul  R2 


Crosshole  and  Downhole  Surveys 


17.  S -wave .  The  results  of  the  crosshole  shear  wave  test  are 
presented  in  figure  9.  A  plot  of  averaged  velocities  for  all  the  tests  versus 
depth  shows  a  steady  increase  in  velocity  with  depth  except  for  an  inversion 
zone  between  8  and  21  ft  (el  0  to  el  -13).  This  data  set  was  divided  into 
seven  velocity  zones  which  adequately  defines  this  trend.  The  interpreted 
layers  and  velocities  are  presented  in  table  2.  The  downhole  S-wave  test 
results  shown  in  figure  10,  agree  well  with  the  crosshole  tests.  For  the 
downhole  test,  the  weakness  of  the  shear  waves  propagating  through  the  deeper 
material  made  the  data  difficult  to  analyze  below  a  depth  of  30  ft  (el  -22). 


Table  2.  Interpreted  Layers  and  S-wave  Velocities 


Crosshole 


Shear  Wave 

Interface  Velocity 
Depth  (el), ft _ 
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Figure  10.  Dovmhole  Shear  Wave  Test  Results 
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18.  P-wave .  The  results  from  the  crosshole  P-wave  test  indicate  a  five 


layer  system,  with  the  interpreted  velocities  and  interfaces  shown  in  table  3 
and  figure  11.  Again,  the  velocities  show  a  general  increasing  trend  with 
increasing  depth.  The  results  of  the  downhole  P-wave  test  are  presented  in 
figure  12.  The  downhole  test  shows  an  increasing  velocity  with  depth  profile 
which  is  also  presented  in  table  3.  The  crosshole,  downhole,  and  surface 
seismic  tests  show  good  agreement  in  the  interpreted  velocities  and  interfaces 
of  the  soil  layers. 


Table  3.  Interpreted  Layers  and  P-wave  Velocities 


Crosshole 

P-Wave 

Interface  Velocity 
Depth  (el) .ft  fps 
1500 

5  (3)  . 

2150 

9  (-1)  . 

5200 

15  (-7)  . 

6200 

20  (-12)  . 

6600 

60  (-52)  . 


Downhole 

P-Wave 

Interface  Velocity 

Depth  (el),  ft _ fps 

1265 

5  (3)  . 

3400 

8  (0)  . 

5000 

15  (-7)  . 

6280 


Interpretation 


19.  To  make  a  meaningful  interpretation  of  the  results  from  the 
different  tests  into  a  generalized  profile  for  the  site,  the  data  were  first 
compared  among  the  different  tests  at  the  same  locations  and  then  the 
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Figure  11.  Crosshole  Compression  Wave  Test  Results 
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data  sets  were  compiled  into  a  composite  profile.  From  these  comparisons  the 
results  were  analyzed  and  the  idealized  P-  and  S-  wave  velocity  profiles  were 
developed.  These  profiles  are  presented  in  figure  13  along  with  the  four 
layer  profile  developed  from  the  soil  samples.  The  first  layer  of  fill  has  a 
P-wave  velocity  range  of  1,500  to  5,200  ft/sec  with  the  5,200  ft/sec  interface 
being  the  top  of  the  water  table  at  a  depth  of  8  ft  (el  0).  This  layer  had  a 
S-wave  velocity  range  of  470  to  600  ft/sec.  The  fill  showed  an  inversion  in 
the  S-wave  velocity  starting  at  the  top  of  the  water  table.  This  inversion 
zone  continues  downward  and  includes  layer  two  which  is  the  beginning  of  the 
original  soil  deposit.  However,  this  layer  shows  an  increase  in  S-wave 
velocity  to  530  ft/sec.  This  inversion  zone  was  given  a  S-wave  velocity  of 
510  ft/sec,  weighted  average  of  the  470  ft/sec  and  530  ft/sec  layers,  for  use 
in  later  soil  modulus  calculations.  Layer  2  did  not  show  an  inversion  in 
P-wave  velocity  which  was  measured  to  be  6,200  ft/sec  because  of  the  material 
being  saturated.  The  third  layer  had  a  P-wave  velocity  of  6,600  ft/sec  and  a 
S-wave  velocity  of  640  ft/sec.  The  fourth  layer  had  a  P-wave  velocity  of 
6,600  ft/sec  and  a  range  of  S-wave  velocities  from  640  to  940  ft/sec.  In 
general,  the  raw  seismic  data  showed  a  some  scatter  which  was  also  evidenced 
in  the  measured  SPT  values.  These  similar  results  suggest  that  the  scatter  in 
the  data  is  mainly  related  to  the  soil  nonhomogenity  and  therefore  the 
velocities  and  properties  reported  are  averages  for  the  soil  at  depths  tested 
with  higher  and  lower  actual  values  existing  depending  on  test  location. 
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.  Composite  Compression  and  Shear  Wave  Site  Profile 


Determination  of  Soil  Material  Parameters 

20.  Determination  of  Shear  Modulus.  Assuming  an  infinite,  homogenous, 
isotropic,  elastic  medium,  and  solving  the  equations  of  wave  motion  for  a 
shear  wave,  which  is  a  wave  confined  to  motion  in  a  plane  in  the  direction  of 
propagation,  results  in  Equation  1  (Richart  et  al ,  1970). 


'/s  (>) 

where 


Vs  -  S-wave  velocity  ft/sec 

2 

G  -  shear  modulus,  lbs/ft 

J>  -  mass  density  of  Soil,  slugs/ft 

o 

6t  -  total  unit  weight  of  soil,  lbs/ft 

2 

g  -  gravitational  constant,  ft/s 


This  equation  states  that  S-wave  velocity  is  dependent  on  the  ratio  of  the 
shear  modulus  to  the  mass  density  of  the  medium.  For  this  particular  site  a 

•5  3 

total  unit  weight  of  100  lb/ft J  was  used  for  unsaturated  soil  and  120  lbs/ft 
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for  saturated  soil.  The  unit  weight  was  converted  to  mass  density;  and,  with 
the  appropriate  shear  wave  velocities  substituted  into  Equation  1,  the  shear 
moduli,  G,  were  computed.  The  results  are  tabulated  in  figure  14. 


21.  Determination  of  Poisson's  ratio.  Again,  assuming  an  elastic  medium, 
if  the  compression-  and  shear-  wave  velocities  are  known,  Poisson's  ratio, 0  , 
can  be  determined  from  the  ratio  of  these  velocities,  Equation  2.  This 
relation  is  given  below  in  Equation  3  (Department  of  the  Army, 1967).  The 
ratios  were  calculated  using  equation  2  and  then  used  in  expression  3  with  the 
results  of  these  calculations  presented  in  figure  14. 


V*  - 

-o  ■- 


Vf/Vf 

2(VS-I) 


(2) 

(3) 


Vp  -  P-wave  velocity,  ft/sec 

Vs  -  S-wave  velocity,  ft/sec 
\)  “  Poiscors'c  ra.iio 


22.  Determination  of  Young's  and  constrained  moduli.  Young's  modulus,  E, 
relates  the  stress  to  the  resulting  strain  when  a  uniform  stress  is  applied  to 
plane  sections  of  a  body  perpendicular  to  the  applied  force  with  the  lateral 
surfaces  free  from  constraint.  The  measurement  of  the  velocity  of  this 
dilational  (P-)  wave  through  the  body  results  in  a  "rod"  velocity  which  with  a 
value  of  the  body's  density,  E  can  be  determined  directly  (Kolsky,  1963). 
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SITE  PARAMETERS 
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Figure  14.  Soil  Moduli  and  Parameters 


However,  with  in  situ  seismic  testing  the  measured  P-wave  is  related  to  the 
constrained  modulus,  M,  and  the  body's  density  as  a  result  of  the  different 
boundary  conditions  that  exist  for  this  type  testing.  The  constrained  modulus 
is  so  named  because  the  lateral  sides  are  constrained  which  more  closely  match 
the  conditions  present  in  an  in  situ  seismic  velocity  measurement.  Knowing 
either  Young's  modulus  or  the  constrained  modulus  the  other  can  be  derived 
using  the  body's  Poisson's  ratio.  In  this  study,  the  shear  modulus  and 
Poisson's  ratio  were  determined,  and  then  the  soil's  Young's  modulus,  E,  was 
calculated  using  Equation  4.  The  constrained  modulus,  M,  was  determined  using 
Equation  5.  These  calculated  modulus  values  are  shown  in  figure  14. 
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Procedure  For  Estimating  Moduli  at  Different  Confining  Stresses 

23.  General .  For  engineering  use  it  is  necessary  to  determine  soil 
moduli  for  the  expected  load  magnitudes  and  loading  conditions  of  the  planned 
structure.  These  moduli  can  be  arrived  at  using  two  approaches  given  that 
certain  in  situ  properties  are  known.  The  first  method  involves  selecting 
properties  at  a  depth  in  which  the  overburden  pressure  is  equal  to  the 
expected  load  of  the  structure  (Department  of  Army,  1967).  This  method  assumes 
that  the  soil  is  homogenous  throughout  the  soil  profile  which  encompasses  the 
depth  of  the  structure  and  the  needed  depth  of  soil  to  equal  the  structure's 
load.  An  alternative  method  allows  an  empirical  constant,  K2 ,  to  be 
determined  for  different  layers  in  the  profile  and  gives  the  advantage  of 
using  soil  parameters  that  can  vary  in  depth  and  therefore  model  the  soil  more 
precisely  at  the  proper  depth. 

24.  Determination  of  K2  Values.  Investigators  have  shown  that  for  sands 
shear  modulus  values  are  strongly  influenced  by  the  confining  pressure,  void 
ratio  or  relative  density,  and  strain  amplitude  (Seed  et  al,  1970).  For  the 
case  of  low  amplitude  strains  <  5n(10"^),  which  is  the  range  of  in  situ 
seismic  wave  velocity  measurements,  a  Gmax  is  determined,  and  the  resulting 
coefficient  in  equation  6,  K2 ,  is  a  K2max,  which  is  also  applicable  to  elastic 
response  loading  conditions.  It  should  be  further  noted  that  these  field 
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velocities  are  determined  from  total  stress  conditions  which  means  that  the 
moduli  are  for  undrained  loading  conditions.  The  expression  that  relates 
Gmax  to  these  factors  is  given  in  Equation  6. 

G  * /Qoofci/lfr'  U) 

lAjnert  __  . 

C~/>7  '  1  -J  ( *  s~i  4  ) 

-  o.az  <r*' 

Gmax  -  low  strain  amplitude  shear  modulus 

K2  -  empirical  proportionality  constant 
relating  Gmax  to  Y  air*' 

Kq  -  coefficient  of  lateral  pressure 
(assigned  a  value  of  0.45) 

2 

Cr “  effective  mean  confining  stress,  lbs/ft 

2 

<TV'  "  effective  vertical  stress,  lb/ft 

It  can  be  seen  that  the  relational  factor  K2  can  be  determined  from  data  at 
hand  as  a  function  of  depth  for  the  soil  profile.  These  resulting  K2's  can  be 
averaged  within  each  layer  with  the  resulting  average  K2  assigned  to  that  zone 
of  material.  This  assigned  K2  can  then  be  used  to  estimate  shear  moduli  at 
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selected  confining  stresses  at  any  point  in  the  profile  by  use  of  Equation  6 
Table  4  given  below  is  derived  from  data  given  in  Seed  et  al .  (1970)  and  can 


be  used  to  estimate  relative  density  for  sands  from  K2  values. 


Table  4. 

Estimating  Relative  Density  of  Sands  from  K2  Values 


Relative  Density  % 

K2 

Loose 

30 

34 

40 

40 

Medium 

45 

43 

60 

52 

Dense 

75 

61 

Very 

90 

70 

Dense 

25.  The  calculated  K2's  for  each  depth  increment  in  the  soil  profile  are 


given  in  figure  14.  A  value  of  Ko-0.45  was  assumed  in  determining  the  mean 
effective  confining  pressure. 


26.  Soil  Moduli  Determinations.  The  use  of  the  K2  parameter  to  estimate 
soil  moduli  for  various  loading  conditions  can  be  illustrated  using  this 
procedure  to  compare  the  moduli  derived  from  the  seismic  and  plate  bearing 
tests.  The  moduli  for  the  expected  load  was  calculated  using  the  sum  of  the 
static  plus  one-half  of  the  cyclic  plate  bearing  loads  (3,140  lbs/ft  )  as  the 
vertical  effective  stress.  The  in  situ  moduli  from  the  seismic  results  were 
obtained  from  the  same  test  elevation  as  the  plate  bearing  test.  These  same 
calculations  were  repeated  to  a  depth  of  30  ft  (el  -22).  A  table  of  moduli 
values  for  this  applied  load  to  a  depth  of  30  ft  and  a  comparison  of  Young's 
modulus  between  the  seismic  and  plate  bearing  tests  are  presented  in 
figure  15.  The  results  from  the  seismic  tests  estimate  a  lower  E  than  the 
plate  bearing  test  results.  The  seismic  results  are  strongly  dependent  on  the 
magnitude  of  the  S-wave  velocity.  The  reported  S-wave  velocity  of  600  ft/sec 
is  a  conservative  average;  however,  if  a  velocity  of  720  ft/sec  is  used,  the 
agreement  is  much  closer.  This  velocity  is  a  raw  average  of  all  S-wave  data 
at  a  depth  of  2.5  ft.  This  result  is  also  presented  in  figure  15  for 


comparison. 


SOIL  MODULI  ESTIMATIONS 


DEPTH 

ft 

LOAD 
lb/  ft2 

K2 

Poisson ' s 
ratio 

G*  E* 

lb/ in'. 

2.5 

3140 

89 

0.410 

27561 

77722 

5.0 

3140 

63 

0.410 

19510 

55017 

7.5 

3140 

51 

0.470 

15793 

46433 

10.0 

3140 

32 

0.490 

9910 

29531 

12.5 

3140 

36 

0.490 

11148 

33222 

15.0 

3140 

34 

0.490 

10529 

31376 

17.5 

3140 

32 

0.490 

9910 

29531 

20.0 

3140 

31 

0.490 

9600 

28608 

22.5 

3140 

46 

0.490 

14245 

42450 

25.0 

3140 

44 

0.490 

13626 

40605 

27.5 

3140 

43 

0.490 

13316 

39682 

30.0 

3140 

41 

0.490 

12697 

37836 

*  Based  on  K2  soil  parameter 


COMPARISON  OF  MODULI  DETERMINATIONS 

Source  Load  K2  Poisson's  G  EM 

lb/ ft*  ratio  lb/  in2. 


Seismic  Test  Procedure 


Vs=600  Vp=1500 , 
3140  89 

ft/sec 

0.410 

27561 

77722 

180678 

Vs=720  Vp=1500 , 
3140  128 

ft/sec 

0.350 

39638 

107023 

171764 

Plate  Bearing  Test 

123000 

Figure  15.  Comparison  of  Moduli  Determinations  from  Seismic 
and  Plate  Bearing  Test  Methods 
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Conclusions 


27.  From  Che  data  that  have  been  presented  the  following  general 
conclusions  are  made. 

a.  The  site  can  be  divided  into  the  following  zones  with  material 
classifications  and  characteristic  P-  and  S-wave  velocities  as  follows: 


Material 

P - wave 

S - wave 

(ft/sec) 

1 

(ft/sec) 

Layer  1:  sand  fill 

i 

SP 

:  15001  -  52002 

470  -  600 

Layer  2:  sand, silty 

i 

1 

| 

SM 

5200  -  6200 

530 

Layer  3:  sand, fine 

■ 

!  sw 

6600 

640  | 

i  Layer  4:  sand, silty 

i 

i 

i 

!  gravely 

1 

i 

i 

SM 

6600 

640  -  940  [ 

I 

1  2 
unsaturated  saturated 

b.  The  SPT  test  results  show  that  the  measured  N- values  have 
considerable  scatter  especially  in  the  upper  10  ft  at  the  site.  The  scatter 
is  probably  due  to  the  naturally  occurring  gravel  and  cobbles  that  were  in 
layer  one  and  buried  construction  materials  in  the  near  surface.  The  average 
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equivalent  (Nl)gQ  blowcount  for  the  fill  is  25  with  the  original  sands  between 
14  and  29  ft  in  depth  (el  -6  to  -21)  having  a  value  of  21.  The  sands  below  a 
depth  of  29  ft  (el  -21)  show  a  decrease  in  the  (N1)^q  blowcount  to  a  value  of 
15.  A  global  site  value  for  the  sands  would  be  22  which  corresponds  to  a 
medium  dense  sand. 

c.  The  Shear  and  Young's  moduli  and  Poisson's  ratio  have  been 
determined  for  each  test  depth  and  tabulated  in  figure  14.  These  parameters 
have  been  grouped  according  to  idealized  layers  and  presented  in  figure  16. 
Throughout  the  entire  profile  the  shear  modulus  ranges  from  7. 7 *(10  )  Ibs/in 
at  the  surface  to  22* (10^)  lbs/in^  at  a  depth  of  60  ft.  Young's  modulus 
varies  from  22* (10  )  lbs/in  at  the  surface  to  66*(10  )  lbs/in  at  the  bottom 
of  the  profile.  The  unsaturated  materials  near  the  surface  have  a  Poisson's 
ratio  of  0.40  with  the  saturated  materials  below  the  water  table  having  an 
expected  value  of  0.49.  There  exists  a  zone  between  10  and  20  ft  in  depth  (el 
-2  to  -12)  where  the  moduli  decrease  due  to  the  inversion  in  the  S-wave 
velocity  profile  at  these  depths. 

d.  A  K2  parameter  was  calculated  for  each  test  depth  and  an  average 
K2  assigned  to  zones  in  the  profile  and  shown  in  figures  14  and  16.  Layer  one 
has  a  K2  of  68  in  the  upper  7  ft  and  33  for  the  remainder.  The  second  layer 
also  has  a  K2  value  of  33.  The  upper  2  feet  of  layer  three  has  a  K2  of  33 
which  increases  to  42  below  that  depth.  The  fourth  layer  ranges  from  a  K2  of 
42  at  a  depth  of  29  ft  (el  -21)  to  a  K2  of  65  at  a  depth  of  60  ft  (el  -52). 
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LAYER  DEPTH  XT  \  SEISMIC  VELOCITY  SOIL  PARMAMETERS 
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Figure  16.  Summary  of  Soil  Moduli  and  Parameters 


e.  The  K2  parameter  was  used  to  estimate  moduli  for  a  load  of  3,140 

2 

lbs/ft  and  then  compared  with  the  plate  bearing  test  results.  The  plate 

2 

bearing  test  estimates  an  E  of  123,000  lbs/in  which  is  higher  than  the  E  of 
o 

78,000  lbs/ft  derived  from  the  K2  parameter. 

f.  Comparing  the  relative  densities  of  the  sands  at  the  site  using 
correlations  based  on  calculated  K2  values  and  (N1)^q  blowcounts  shows  good 
agreement  except  that  the  K2  derived  relative  density  estimates  are  lower  for 
the  sands  between  the  10  and  30  ft  depth  range  (el  -2  to  -  22).  These 
relative  density  estimates  are  shown  in  figure  16. 
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APPENDIX  A 
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APPENDIX  B 


PLATE  BEARING  TEST 
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